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To execute a trade, participants in electronic equity markets may choose to submit limit orders or market 
orders across various exchanges where a stock is traded. This decision is influenced by the characteristics 
of the order flow and queue sizes in each limit order book, as well as the structure of transaction fees and 
rebates across exchanges. We propose a quantitative framework for studying this order placement problem by 
formulating it as a convex optimization problem. This formulation allows to study how the interplay between 
the state of order books, the fee structure, order flow properties and preferences of a trader determine the 
optimal placement decision. In the case of a single exchange, we derive an explicit solution for the optimal 
split between limit and market orders. For the general problem of order placement across multiple exchanges, 
we propose a stochastic algorithm for computing the optimal policy and study the sensitivity of the solution 
to various parameters using a numerical implementation of the algorithm. 
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1. Introduction 



In todays' automated, electronic financial markets, the trading process is divided into several 
stages, each taking place on a different time horizon: portfolio allocation decisions are usually 
made on a monthly or daily basis and translate into trades that are executed over time intervals 



of several minutes to several days. Existing studies on optimal trade execution (Bertsimas and Lo 



1998 Almgren and Chriss||2000 ) have investigated how the execution cost of a large trade may be 



reduced by splitting it into multiple orders spread in time. Once this order scheduling decision is 
taken, one still needs to specify how each individual order should be placed: this order placement 
decision involves the choice of an order type (limit order, market order), order size and destination, 
when multiple trading venues are available. Orders are filled over short time intervals of a few 
milliseconds to several minutes and the mechanism through which orders are filled in the limit 
order book are relevant for such order placement decisions. Market participants need to make such 
decisions thousands of times each day, and their outcomes have a large impact on each participant's 
transaction cost as well as on aggregate market dynamics. 



Early work on optimal trade execution (Bertsimas and Lo 1998, Almgren and Chriss 2000) did 
not explicitly model the process whereby each order is filled, but more recent formulations have 



tried to incorporate some elements in this direction. In one stream of literature (see Obizhaeva and 



Wang (2005), Alfonsi et al. (2010), Predoiu et al. (2011)) a trader is restricted to using market 



orders whose execution costs are given by an idealized order book shape function. Another approach 



is to model the process through which an order is filled as a dynamic random process (Cont (2011 ), 



Cont and De Larrard (2011)) and thus formulate the optimal execution problem as a stochastic 



control problem: this formulation has been studied in various setting with limit orders (Bayraktar 



and Ludkovski (2011 ), Gueant and Lehalle (2012)) or limit and market orders (Guilbaud and Pham 



2012, Huitema 2012) but its complexity makes it intractable unless restrictive assumptions are 



made on price and order book dynamics. 

In the present work, we adopt a simpler, more tractable approach: assuming that the trade 
execution schedule has been specified, we focus on the task of filling each order. Decoupling the 
scheduling problem from the order placement problem leads to a more tractable approach which 
is closer to market practice and allows us to incorporate some realistic features which matter for 
order placement decisions, while conserving analytical tractability. 

Individual order placement and order routing decisions play an important role in modern financial 
markets. Brokers are commonly obliged by law to deliver the best execution quality to their clients 
and empirical evidence confirms that a large percentage of market orders in the U.S. and Europe 



is sent to trading venues providing lower execution costs or smaller delays ( Boehmer and Jennings 
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2007 Foucault and Menkveld||2008 ) . Market orders gravitate towards exchanges with larger posted 



quote sizes and low fees, while limit orders are submitted to exchanges with high rebates and lower 
execution waiting times (see Moallemi et al. ( 2011| )). These studies demonstrate how investors' 
aggregate order routing decisions have a significant influence on market dynamics, but a systematic 
study of the order routing problem from the investor's perspective is lacking. A reduced-form model 



for routing an infinitesimal limit order to a single destination is used by Moallemi et al. (2011), 



while Ganchev et al. (2010) and Laruelle et al. (2009) propose numerical algorithms to optimize 
order executions across multiple dark pools, where supply/demand is unobserved. To the best of 
our knowledge this paper is the first to provide a detailed treatment of investor's order placement 
decision in a multi-exchange market, unified with the market/limit order choice. 

Our key contribution is a quantitative formulation of the order placement problem which takes 
into account multiple important factors - the size of an order to be executed, lengths of order 
queues across exchanges, statistical properties of order flows in these exchanges, trader's execution 
preferences, and the stucture of liquidity rebates across trading venues. Our problem formulation is 
tractable, intuitive and blends the aforementioned factors into an optimal allocation of limit orders 
and market orders across available trading venues. Order routing heuristics employed in practice 
commonly depend on past order fill rates at each exchange and are inherently backward-looking. In 
contrast, our approach is forward-looking - the optimal order allocation depends on current queue 
sizes and distributions of future trading volumes across exchanges. When only a single exchange is 
available for execution, this order placement problem reduces to the problem of choosing an optimal 
split between market orders and limit orders. We derive an explicit solution for this problem and 
analyze its sensitivity to the order size, the trader's urgency for filling the order and other factors. 
Similar results are also established in a case of two trading venues under some assumptions on 
order fiow distributions. Finally, we propose a stochastic approximation method for solving the 
order placement problem in the general case and demonstrate its efficiency through examples. Our 
numerical examples demonstrate that the use of our optimal order placement method allows to 
substantially decrease trading costs in comparison with various 'naive' order placement strategies. 

An important aspect of our framework is to account for execution risk, through the incorporation 
of a penalty for under- or over-filling an order. This penalty is high for time-sensitive executions 
or when it is costly to catch up on the unfilled portion of the order. Although market orders are 
executed at a less favorable price, it becomes optimal to use them when execution risk is a primary 
concern. Optimal limit order sizes are strongly influenced by total quantities of orders queueing for 
execution at each exchange and by distributions of order outflows from these queues. For example, 
if at one of the exchanges the queue size is much smaller than the expected future order outflow, 
it is optimal to place a larger limit order there. Finally, the total order size to be filled plays an 
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important role - limit orders are used predominantly to execute small order sizes and market orders 
are used for medium and large orders. The amount that can be realistically filled with a limit 
order at each exchange is naturally constrained by the corresponding queue size and order outflow 
distribution, so the share of market orders in the optimal allocation increases as the total order 
size increases. We find that the optimal order allocation almost always splits the total quantity 
among all available exchanges, suggesting that there is a benefit in having multiple markets. 

Section [2] describes our formulation of the order placement problem and shows that it has a 
global optimum. In Section [3] we derive an optimal split between market and limit orders for a 
single exchange. Section |4] analyzes the general case of order placement on multiple trading venues. 
Section [5] presents a numerical algorithm for solving the order placement problem in a general case 
and our simulation results, and Section [6] concludes. All proofs are presented in the Appendix. 

2. The order placement problem 

Consider a trader who has a mandate to buy S shares of a stock within a (short) time interval 
[0,T]. The deadline T may be a fixed horizon (e.g. 1 minute) or a stopping time (triggered by 
market activity) . To gain queue priority the trader may immediately submit K limit orders of sizes 
Lk to various exchanges k = 1, . . . , K or submit one market order of size M. The trader's order 
placement decision is thus summarized by a vector X = (M, Li, . . . , Lk) G whose components 

are non- negative (only buy orders are allowed). Our objective is to define a meaningful framework 
in which the trader may choose the various possibilities for this order placement decision. 

We focus on limit order placement and execution and assume that a market order of size M can 
be filled immediately and with certaintjj^ Limit orders with quantities (Li, . . . , Lk) join queues of 
(Qi, . . . ,Qk) pre-existing limit orders at the best bids of K exchanges, where Qk > 0. To simplify 
the notation, we make an assumption that all K available bid queues are lined up at the best bid 
price, but it is easily relaxed. Denote by {x)+ = max(a;, 0). If Lk is constant within [0,T], the amount 
purchased with a limit order on exchange k by time T is equal to {^,k — Qk)+ — (Cfe ~ Qk — Lk)+, 
where £,k = Ck + is an order outflow from the front of /c-th bid queue, consisting of G [0, Qk\ 
cancelations of pre-existing orders from that queue and trades with contra-side marketable 
orders reaching that queue. We specifically note that limit order fill amounts are random, and we 
allow for partial fills. The total amount A{X,^) bought by the trader by time T with all of his 
orders is a function of the order allocation X and an overall bid queue outflow ^ = (Ci, • • • , C-R")- 

^ This assumption is reasonable if S is small relative to the prevailing market depth. Under the assumption of 
immediate and certain market order execution it is easy to show that sending market orders to exchanges with high 
fees is always sub-optimal. We therefore consider a single exchange (with the smallest liquidity fee) for the purpose 
of sending a single market order. 
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K 

A{X, ^) = M + Y, ((6 - Qk)+ - (6 -Qk- Lk)+) (1) 

The total price of this purchase is divided into a benchmark cost paid regardless of trader's 
decisions, computed using a mid-quote price level, and an execution cost given by 

K 

{s + f)M - Y,{s + nM, - Qk)+ - [ik -Qk- Lk)+), (2) 
fe=l 

where s is a half of the bid-ask spread at time 0, / is the lowest available liquidity fee and rk,k = 
1,...,K are liquidity rebates for all exchanges. The trader can reduce the execution cost by sending 
more limit orders, but this leads to a risk of falling behind the target quantity S because their 
fills are random. To capture this execution risk we include, in the objective function, a penalty for 
violations of target quantity in both directions: 

A„ {S - A{X, 0)+ + K {A{X, (3) 

where A„,Ao are marginal penalties for, respectively falling behind or exceeding the execution 
target S. These penalties are motivated by a correlation that exists between limit order executions 
and price movements (so-called adverse selection). If A{X,^) < S, the trader has to purchase 
the remaining S — A{X,$,) shares at time T with market orders. Adverse selection implies that 
conditionally on the event {A{X,$,) < S} prices have likely moved up and the transaction cost 
of market orders at time T is higher than their cost at time 0, i.e. A„ > s -|- /. Alternatively, if 
A{X, S the trader experiences buyer's remorse - conditionally on this event prices have likely 
moved down and he could have achieved a better execution by being more patient. Besides adverse 
selection, parameters A„,Ao may reflect trader's execution preferences. For example a trader with 
a positive forecast of short-term returns may prefer to trade early with a market order and set a 
larger value for A„. 

Problem 1 (Optimal order placement problem) An optimal order placement is a vector 
X* e M++^ solution of 

mm E[v{X,0] (4) 

where 

K 

v{X,i):={s + f)M-Y,{s+rkmk-Qk)+-{^k-Qk-Lu)+) + Xu{S-A{X,£)))^ + X,{A{X,i)-S)^ 
fe=i 

(5) 

is the sum of the execution cost and penalty for execution risk. 

We will denote V{X) = E[v{X,^)]. We begin by assuming certain economically reasonable restric- 
tions on parameter values. 
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Assumptions 

Al A„ > 0, Ao > 0: the trader is penalized for falling behind or exceeding the target quantity 

A2 Ao > s + max{rfc} and Ao > — (s + /): it is suboptimal to exceed the target quantity S regardless 

k 

of fees and rebates 

A3 minjrfc} + s > 0: even if some are negative, limit orders still reduce the execution cost. 

k 

Proposition [T] below shows that it is not optimal to submit limit or market orders that are a priori 
too large or too small (larger than the target size S or whose sum is less than S). Proposition [2] 
guarantees the existence of an optimal solution. 

Proposition 1 Consider C - a compact convex subset o/M+^^ defined by 



C = |xeM++^ 0<M<S, 0<Lk<S-M,k = l,...,K, M + ^ 
I fc=i 



Under assumptions A1-A3 for any X ^ C, 3X' G C with V{X') < V{X). Moreover, if 
min{P(^fc > Qk + S)} > 0, the inequality is strict: V{X') < V{X). 

k 

The penalty function ([s]) implements a soft constraint for order sizes and effectively focuses the 
search for an optimal order allocation to the set C. Specific economic or operational considerations 
could also motivate hard constraints, e.g. M = or XlfLi -^^ ~ ^- Such constraints can be easily 
included in our framework but absent the aforementioned considerations we do not impose them 
here. 

Proposition 2 Under assumptions A1-A3, V{X) is a convex function on M;^'*'^, it is bounded 
below and has a global minimizer X* £ C. 
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3. Choice of order type: limit orders vs market orders 

To highlight the tradeoff between hmit and market order executions in our optimization setup, we 
first consider a case when the asset is traded on a single exchange, and the trader has to choose an 
optimal split between limit and market orders. Since K = 1, we suppress a subscript 1 throughout 
this section. 

Proposition 3 (Single exchange: optimal split between limit and market orders) 

Assume that ^ has a continuous distribution and (Al-AS) hold. Denote 

, A 2g + / + r —A 2s + / + r 

If < Xu, the optimal allocation is {M*,L*) = (0,5). // A„ > A„, the optimal allocation is 
(M^jL*) = (S*, 0). If Xu & (A„,A„), the optimal allocation is: 



Xu + s + r ^ 

where F{-) is a cumulative distribution function of the bid queue outflow ^. 

In the case of a single exchange. Proposition [T] implies that M* + L* = S, therefore there is no 
risk of exceeding the target size and Aq does not affect the optimal solution. The trader is only 
concerned with the risk of falling behind the target quantity, and balances this risk with the fee, 
rebate and other market information. The parameter A„ can be interpreted as trader's urgency to 
fill the orders, and higher values of A„ lead to smaller limit order sizes, as illustrated on Figure [ij 
In contrast, the optimal market order size increases with A„. 

The optimal split between market and limit orders depends on the ratio x^^^g^^^-^ which balances 
marginal costs and savings from a market order. It also depends on the distribution F and the 
queue length Q - keeping all else constant, a trader would submit a larger limit order if its execution 
is more likely and vice versa. The optimal limit order size decreases with A„ as it becomes more 
expensive to underfulfill the order and increases with / as market orders become more expensive. 
Another interesting feature is that L* is fully determined hy Q, F and pricing parameters s, r, /, A„, 
while M* increases with S. As a consequence of this solution feature, as the order size S increases, 
a larger fraction ^ of that order is executed with a market order. The solution (M*,L*) depends 
on the entire distribution of ^ and not just on its mean, as illustrated on Figure [l] for a pair of 
exponential and Pareto distributions with equal means. 
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4. Optimal routing of limit orders across multiple exchanges 

When multiple trading venues are available, dividing the target quantity among them provides 
better execution quality by reducing the risk of not filling the order. However, sending too many 
orders leads to an undesireable possibility of exceeding the target size. Proposition|4]gives a criterion 
for optimality of an order allocation X* = {M* ,L\, . . . ,L*i^) that balances these risks. 

Proposition 4 Assume (A1-A3), also assume that the distribution of ^ is continuous, 

( 2s + f + 1 

max{Ffc(Qfe + S')} < 1 and A„< max <^ — {s + rk)\. Then: 

^ ^ {_ rk[Qk) J 

1. Any optimal allocation X* has a positive market order quantity M* > if 

2s + f + maxjrfc} 

A„ > — y ^ - (s + m_ax{rfe}) (7) 

< Qk} 



2. Any optimal allocation X* has a positive limit order quantity L* > if 



r\{Ck<Qk} 



3. // ([7])-([8]) hold for all exchanges j = 1, . . . , K , a necessary and sufficient condition for opti- 
mality of an order allocation X* is that it solves the following equations: 



K 



M* + Y^ (fe - Qk)+ - iCk -Qk- Ll)+) < S 
\ fe=i 

^ K 

M* + Y^ (fe - Qk)+ - (6 -Qk- Ll)+) < S 



S + f + Xo 
Xu + 



(9) 



k=l 



ij>Q,+L* 



_ Aq - (s + rj- 

Ati + Ao 
3 = 1,. ..,K 



(10) 



Equations ([9])-(10) show that an order allocation is optimal as long as it sets the probabilities of 
falling behind the target quantity equal to specific thresholds computed with pricing parameters. 
When the number of exchanges K is large, the probabilities in Q-flO) are difficult to compute in 
closed- form. However, before turning to numerical procedures we investigate how these equations 
can be solved in a tractable case of two exchanges. 
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0.025 0.035 0.045 0.055 0.065 0.075 



Figure 1 Optimal limit order size L* for one exchange. The parameters for this figure are: Q = 2000, S — f 000, s = 
0.02, r — 0.002,/ = 0.003. Colors correspond to different order outflow distributions - exponential with 
means 2200 and 2500 and Pareto with mean 2200 and a tail index 5. 



Corollary Consider the case of two exchanges with £,i,£,2 that are independent and have continuous 
distributions. If 



1. max{Ffe(Qfe + 5)}<l, 



2. A„ < max 

k=l,2 



2s + f + rk 



3. Fi(Qi)<l-^^,F2(g2)<l 



s + ri 



2s + f + maxjrfc} 



Ao Ao 
then there exists an optimal order allocation X* = (M*,L^,L2) G int{C} and it verifies 



Au + Ao 
-1 / Ao - (s + ra) 



Ll = Q, + S-Ar -F- 
Fi(gi + Ll)F^{Q., + S-Ar- LI) + 



(11a) 
(lib) 



F2{Qi +Q2 + S-AF- x,)dF,{x,) = ^\ + (11c) 



Ql + S-M* -L* 

where Fi{-),F2{-) are the cd/ 0/^1,^2 respectively. 

In this solution optimal limit order quantities L*,L2 are linear functions of an optimal market 
order quantity M*. When (lla)-(llb) are substituted into ( |llc ) we obtain a (non-linear) equation 
for M* , which can be solved for a given distribution of (^1,^2)- 
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Example If ^1,^2 are exponentially distributed with means //i,//2 respectively, then an optimal 
order allocation is given by: 

f M* = Qi + Q2 + 5 - z 

r* ^ , 1 fK + s + ri 

L^ = z-Qi+H2 log 

Ll = z- Q2+ fii log 
where z is a solution of a transcedental equation: 



Xu + Xo J (12) 
A„ + s + r2 ' 

Xu + Xo 



1 , 1 ,'(A„ + s + ri)(A„ + s + r2)\ z A„-(s + /) z 

1 + log ^^ — r^r~^9 +~ = — — — e^^ , if ^1 =/i2 = ^ (13) 



or 



{Xu + XoY J 11 A„ + Ao 

^1 -^/A„ + s + ri\ ''I ^2 -^/A„ + s + r2\ "2 A„ - (s + /) 

-e ''I — — — — H e "2 — — — — = — - — — , if ^1 / ^l2 



\ A„ + Ao / IJ'2 ~ fJ-i V An + Ao / Xu + Ao 

(14) 

Similarly to the case of one exchange, in this example an optimal market order size M* is an 
increasing linear function of queue sizes Qi,Q2 and the target quantity S, while optimal limit 
order sizes L* are decreasing functions of the corresponding queue sizes Qi. As in the case of one 
exchange, the optimal limit order sizes do not depend on the target quantity, but the optimal 
market order size increases with it. In addition we note that each L* depends on the order flow 
distribution on both exchanges through /ii 2 and z. 

5. Numerical solution of the optimization problem 

Computing the objective function in the order placement problem Q or its gradient at any point 
requires calculating an expectation (a multidimensional integral) which, aside from special exam- 
ples, is generally not analytically tractable. Stochastic approximation methods, developed specif- 
ically for problems where the objective function is an expectation, turn out to be very useful for 
this problem. We propose in this section a procedure for computing the optimal order placement 
policy which does not require specifying an order outflow distribution. 

We begin by briefly discussing the stochastic approximation approach and the specific method 
used here. Consider an objective function V{X) = E,[v{X,^)] to be minimized and denote by 



g{X,(^) = VvlXjS^) where the gradient is taken with respect to X. The Robbins and Monro (1951) 



stochastic approximation algorithm tackles the problem in the following way: 
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1: Choose Xq and a sequence of 'step sizes' {7^}; 
2: for n = l,...,N do 

3: Simulate an independent random variable ^" with distribution F 

4: Set X„ = X„_i -7„ff(X„_i,e) 
5: end 

This algorithm produces an estimate X* = X^ , which converges to the optimal point X* as — t- 00 
under some technical assumption^ The sequence of constants {7„} affects the rate of convergence. 
This sensitivity can be overcome by using for example the robust stochastic approximation of 



Nemirovski et al. ( 2009 ) which follows the same steps as above with a constant step size 7,1 = 7 



and uses an average of iterates X* = ^n=i -^n instead of X[^ as an estimate of the optimal point. 
Under some weak assumptions this method achieves a performance bound V{X*) — V{X*) < 



for a finite A^, where D = max \\X — X'\\2 and M = y'max E with || • || denoting L2 

nornj^ Multiplying an optimal step size 7 by a constant 9 >0 leads to a performance bound of the 
same order of magnitude max{0, 0"^}^^, i.e. the method is "robust" to step size misspecifications. 



For our problem we can further bound 



K 



D<VkS and M<i{s + f + X^ + X,f + ^{s + rk + Xu + Xo) 



k=l 




We assume that on each iteration X„ G int{C} - this is enforced by rescaling X„ when needed and 
does not affect the convergence - then the stochastic gradient in our problem is given by: 



/ s + / — A„l{A(js:„,e)<s} + Aol{A{x„,5)>5} \ 



{s + ''l)l{^i>Qi+£i „} '^"l{A(X„,e)<S,ei>Qi+Li,„} + {A(X„,0>S,ii>Qi+Li,„} 
V ^^^^'^'^^^{iK>QK+LK,^} ^ •^"-'^{.4{X„,C)<S,5,f>Q;f + Lj^.„} +'^°^{A(X„,0>S,«K>QK+i/f.„}/ 

Note that (jr(X„,^) depends on random variables ^ only through indicator functions, which have 
simple economic meaning. For example 1{a(x„,^)<5} = 1 if the target quantity was not fully executed 
on the n-th step and '^^^^.>Qf.+L^. | = 1 if there was an opportunity to execute a larger limit order 
at exchange k on that step. This leads to a simple interpretation of numerical iterations - at each 
step order sizes are increased or decreased depending on whether or not the executed quantity 



2 



Specifically, if V{X) need to be strongly convex, twice continuously differentiable and X* £ int{X}, where A" is a 



nonempty bounded closed convex set, then E[V(X„) — — > 0. We refer to Kushner and Yin (20031 for more 
details. 

^ The method assumes that min{V^(X)} is sought, where V{X) is a well-defined and finite-valued expectation for 

every X £ X and A" is a non-empty bounded closed convex set. Moreover V{X) needs to be contirmous and convex 
on X. The optimal step size is 7= ^^^^ . 
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was smaller or larger than the target size and whether or not a larger limit order could be filled 
at a given exchange. If a model for is available, one can use it to simulate and compute a 
numerical solution that takes into account specific order fiow assumptions (e.g. forecasts of future 
trading volume). Alternatively, one can use past order fill data to compute indicator functions in 
g{Xn,0 obtain a non-parametric numerical solution for the order placement problem (using 
the empirical distribution of past order fills instead of assuming a functional form for F). 

We analyze the numerical stability and convergence of estimates X* by comparing them with an 
analytical solution in the case of one exchange. For this computation we use Q = 2000 shares, ^ ~ 
Pois{^T), fj, = 2200 shares per minute, T = 1 minute and S = 1000 shares. The pricing parameters 
(in dollars per share) are s = 0.02, r = 0.002,/ = 0.003 and fall in a typical value range for US 
equities. Finally, the penalty costs (also in dollars per share) are set to Ao = 0.024, A„ = 0.026. 
According to Q the optimal allocation (M*,L*) = (728,272) shares. Numerical estimates X* 
were computed for five initial points Xq and different number of iterates N in the stochastic 
approximation, using a step size 7 = ^/KS I N{s + f + Xu + Ao)^ + iV ^ (s + + A„ + Aq)^ I 
For each choice of Xq and N we simulated an additional L = 1000 observations of f to estimate 

1 L n n 

the objective values V{X) with sample averages W{X) = — ^t;(X, ^j)). Figures 2 and 3 show 

^4=1 U U 

that estimates converge to X* regardless of Xq. When Xq = X*, estimates remain close to that 

point. Convergence is also quite fast - after as few as 50 samples the algorithm can be within 2% 

of the optimal objective value. In the worst case of initial points on the boundary it can take a few 

thousand samples to converge. It is also worth noting that convergence in terms of the objective 

value occurs significantly faster than convergence in terms of the order allocation vector. 

We also estimated savings from optimal order placement and from dividing a limit order among 

multiple exchanges. Denote a pure market order allocation by Xm = (5, 0, . . . , 0), a single limit order 

S S S 

allocation by X^ = (0, S*, 0, . . . , 0) and an equal split allocation by Xe = {— , — , . . . , — ). 

K -\- 1 K -\- 1 K -\- 1 

Table 1 presents outputs from the numerical algorithm with Xq = Xe,N = 1000, L = 1000 for 

different order sizes S and number of exchanges K = 1, . . . ,5. The parameters s, f, r, A„, Ao are same 

as in the previous simulation and exchanges are identical replicas of each other: r^. = r, = Q and 

Cn,fc ~ Pois{iiT) are i.i.d. copies of each other, where k = 1, . . . , K, n = l,. . . ,N. Order allocations 

produced by stochastic approximation clearly outperform the naive benchmarks, especially when 

a target quantity S is relatively small and cost savings of limit orders can be fully captured. 

Comparing W{Xl) and W{Xe) we also see that splitting a limit order across multiple exchanges 

can be very advantageous when limit order fills are independent. Since multiple exchanges in this 

example are copies of each other, the algorithm splits the total limit order amount equally among 

S 

them. This is not the same as the allocation Xe because the latter sets a market order size to — , 

K + 1 
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which may be too large or too small depending on S and the number of available exchanges. Another 
interesting feature of the numerical solution X* is a tendency to oversize the total quantity of limit 
orders, which is clearly observed for S = 1000, 5000 and K = 4,5. This may be a consequence of 
assumed independence between - by submitting large orders to multiple exchanges the algorithm 
reduces the probability of falling behind the target quanitity with a relatively low probability of 
exceeding it. 




9^4 

10 10 10 

Num. samples 



Figure 2 Convergence of objective values to an optimal point for different inital points. 
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K 


(M* 


LI 


LI 


^3 


LI Ll)/S 


W{Xm) 


W{Xl) 


W{Xe) 


W{X*) 


5 = 500 


1 


0.481 


0.519 








11.50 


3.36 


2.79 


2.76 


2 


0.034 


0.601 


0.615 






11.50 


3.48 


-2.86 


-6.00 


3 


0.003 


0.438 


0.433 


0.421 




11.50 


3.39 


-5.22 


-10.56 


4 


0.002 


0.280 


0.277 


0.273 


0.264 


11.50 


3.52 


-6.44 


-10.84 


5 


0.001 


0.198 


0.215 


0.224 


0.206 0.214 


11.50 


3.35 


-7.24 


-10.91 


5" = 1000 


1 


0.713 


0.287 








23.00 


16.30 


14.80 


14.20 


2 


0.484 


0.343 


0.338 






23.00 


16.43 


5.88 


5.48 


3 


0.268 


0.338 


0.334 


0.336 




23.00 


16.29 


-3.16 


-3.61 


4 


0.055 


0.316 


0.313 


0.351 


0.333 


23.00 


16.48 


-9.27 


-12.12 


5 


0.003 


0.309 


0.300 


0.300 


0.309 0.321 


23.00 


16.44 


-12.88 


-19.76 


5" = 5000 


1 


0.839 


0.161 








115.00 


120.33 


112.83 


107.76 


2 


0.747 


0.192 


0.192 






115.00 


120.44 


105.86 


99.65 


3 


0.693 


0.189 


0.189 


0.189 




115.00 


120.40 


97.35 


90.71 


4 


0.650 


0.186 


0.186 


0.186 


0.186 


115.00 


120.37 


88.64 


81.89 


5 


0.614 


0.167 


0.167 


0.167 


0.167 0.167 


115.00 


120.40 


79.63 


72.93 



Table 1 Savings from order splitting. 



To illustrate the structure of a numerical solution we performed a sensitivity analysis with K = 2 
exchanges and parameters Qi = Q2 = 2000, S = 1000, ^1,2 ~ Pois{fii^2T), Hi = 2600, /is = 2200, T = 
1, s = 0.02, ri = r2 = 0.002, / = 0.003, A„ = 0.26 and Aq = 0.24. Varying some of these parameters 
one at a time we plot the numerical solution X* after = 1000 iterations, together with an 
analytical solution for a single exchange. The results are presented on Figures |4] and [5] Similarly 
to the single exchange case, limit order sizes on two exchanges Li,L2 decrease and market order 
size M increases as the penalty A,j increases. Increasing the half-spread s, the rebate ri or the fee 
/ makes a limit order on exchange number one more attractive, so Li increases and M decreases. 
Because the penalty A^ is large in this example, execution risk is more important than fees and 
rebates, therefore the queue size Qi and the order outflow mean have a much stronger effect 
on the optimal solution than ri. Both decreasing the Qi and increasing /ii make a limit order fill 
more likely at exchange number one and Li increase fj5 Finally, as in the case of a single exchange, 
the target size S has a strong effect on the optimal order allocation. Only limit orders are used 
while S is small, but as it becomes larger it is difficult to fill that amount solely with limit orders 
and the optimal market order size begins to grow to limit the execution risk. 



* The observed drop in Li for large fj,i and small Qi is a feature of this example, we were not able to replicate it for 
other distributions of ^. 
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Figure 5 Sensitivity analysis for a numerical solution X* — {M,Li,L2) with two exchanges and an optimal 
solution (M",L") with the first exchange only. 
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6. Conclusion 

We formulated and solved the problem of placing a small batch of orders on multiple trading 
venues. In the case when there is a single exchange we derived an optimal split between a limit and 
a market order sizes. For more general cases, we proposed and tested a stochastic approximation 
algorithm that is shown to quickly converge to an optimal point. We explored the properties of an 
optimal order allocation policy and showed that splitting an order across multiple exchanges can 
lead to a substantial reduction in transaction costs. 



Appendix. Proofs 

Proposition 1 Consider C - a compact convex subset o/M+"'"^ defined by 



K 



C = lXeRl+^ 0<M<S, 0<Lk<S-M,k = l,...,K, M + ^Lfe>5 



k=l 



Under assumptions A1-A3 for any X ^ C, 3X' G C with V{X') < V{X). Moreover, if 
min{P(^fe >Qk + S)} > 0, the inequality is strict: V{X') < V{X). 

k 

Proof: First, for any allocation X that has M > S, we automatically have A{X) > S and we can 
show that the (random) cost and penalty of X is larger than those of X""'^'^ = (5, 0, . . . , 0) G C: 



K 



v{X, - v{X",0 = {s + f){M-S)-J2i^ + r^m, - Qfe)+ - -Qk- Lk)+)+ 

k=l 

x„(^M-s+J2 (fe - Qk)+ - (6 -Qk- Lk)+)^ = 

K 

{Xo + s + f){M -S) + ^(A„ - s - r^mk - Qk)+ - (6 -Qk- Lk)+) > 0, 



fe=i 



which holds for all random ^. Therefore, V{X) > y (^X""'^'"'') . Similarly, for any allocation X with 
Lk>S-M define a diflFerent allocation X' by M' = M, L'- = Lj^j =^ k and L'^ = S- M. Then 
v{X, ^) — v{X',^) = on the event B = {oo\^k{^) <Qk + S — M}. On its complementary event 

v{X, - v{X',0 = -(« + rkMk -Qk-S + M)+- (^fe -Qk- 
+Ao((6 -Qk-S + M)+- {ik -Qk- Lk)+). 

Therefore 



+ E 



(Ao -{s + rfe))((6 -Qk-S + M)+- (C, -Qk- Lk)+)\B'^ HB") > 
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with a strict inequality if P(i?^) > 0. If X' ^ C, we can continue truncating limit order sizes L'- > 
S — M' following the same argument. Each time the truncation does not increase the objective 
function and finally we obtain X" £ C, such that V{X") < V{X). 

Next, if X is such that M - Xlfcli < S define s = S - M - Xlf^i Lk, take M' = M,L'^ = 
Lk,k = 1, . . . , K — 1 and L'j^ = + s. Then, on the event B = ^uj\S,k{uj) < Qk + we have 
v{X,^) = t;(X',^). However, on the event B'^, 



v{X,0-v{X\0 = {s + rKmK-QK-lK)+-{^k-Qk-lK-s) + ) + \u{{U-QK-LK)+-{ik-Qk-LK-s] 



therefore 



V{X)-V{X')=¥. v{X,i)-v{X',i)\B P(5)+E v{X - v{X' ,i)\B^ ¥{B^) 



+ E 



(A„ + {s + rk)mK -Qk- Lk)+ - (6 -Qk- Lk - s)+)\B' P(S=) > 



with a strict inequality if P(i?^) > 0. 



Proposition 2 Under assumptions Al-AS, V(X) is a convex function on M^^^, bounded from 
below and admits a global minimizer X* G C. 

Proof: First, note that {^k — Qk)+ — {ik — Qk — Lk)+ are concave functions of Lfe. Therefore, A{X,^) 
is concave as a sum of concave functions. Similarly, the cost term in v{X,^) is a sum of convex 
functions, as long as r^ > — s, k = 1, . . . ,K and is itself a convex function. Second, since 5* — A{X, ^) 
is a convex functon of X, and the function h(x) = A„(x)+ — Ao(— x)+ is convex in x for positive 
A„,Ao, so the penalty term h {S — A{X , is also convex. 

If Ao > s + maxjrfc} the function V{X) is also bounded from below since v{X,S,) > —{s + 

k 

max{rk})S. 

k 

Finally, since V{X) is convex, it is also continous and reaches a local minimum Vmin on the 
compact set C at some point X* G C. By convexity, K«m is a global minimum of V{X) on C. 
Moreover, since Ao > s + max{rfe}. Proposition jlj guarantees that Vmin < ViX) for any X ^C, so 
Vjnin is also a global minimum of V{X) on M^^^. ■ 



18 



Cont and Kukanov: Optimal order placement in limit order markets 



Proposition 3 Consider the case of a single exchange where ^ has a continuous distribution. 

(s + r). 



n , h \ A 2s + f + r 2s + / + r 

Denote by A,j = ^ — [s + r) and A„ 



F{Q + S) 
Then, under Assumptions (A1-A3): 



HQ) 



• If ^ the optimal allocation is {M*,L*) = (0,5). // A„ > A„, the optimal allocation is 
{M\L*) = {S,0). 

• If Xu& (Au,A„), the optimal allocation is given by 

Proof: By Proposition [l] there exists an optimal split (M*,L*) gC between limit and market 
orders. Moreover for K = 1 the set C reduces to a line M* + L* = S* so it is sufficient to find M*. 
Restricting L = S-M implies that {A{X,^) > S} = 0, {A{X,(,) < S,C> Q + L} = 0, and we can 
rewrite the objective function as 



ViM) =E[is + f)M-is + r)i{C-Q)+-i^-Q-S + M)+) + 

A„(5-M-((e-Q)+-(e-Q-5 + M)+))))^ 



(15) 



For M € {0,S) the expression under the expectation in (15) is bounded for all ^ and differentiable 
with respect to M for almost all i^, so we can compute V'{M) = '''^Jj^^ by interchanging the order 



of differentiation and integration (see e.g. Aliprantis and Burkinshaw (1998), Theorem 24.5): 



V\M)=E 



Note that if A„ < 



s + f+{s + r)l 

2s + f + r 



Xul{i<Q+s-M}\ =2s + f + r-{s + r + Xu)F{Q + S - M) 

(16) 

{s + r), then V'{M) > for M G (0,5) and therefore V is 

F[Q + S) 

non-decreasing at these points. Checking that V{S) — V{0) > (s + / — A„)5 + (A^ + s + r)5(l — 

2s + f + r 

F{Q + 5)) > we conclude that M* = 0. Similarly, if A„ > — {s + r), then v{M) < 

I kQ) 

for all M G (0, S) and V{M) is non-increasing at these points. Checking that V{S) — V{^) < 
{s + f - Xu)S + (A„ + S + r)S{l - F{Q)) < we conclude that M* = S. Finally, if A„ is between 
these two values, 3e > 0, such that V'{e) < 0,V'{S — e) > and by continuity of V there is a 
point where V'{M*) = 0. This M* is optimal by convexity of V{M) and ^ solves equations 
v{M*,C) = 0,L* = S - M*. U 



Proposition 4 Assume (A1-A3), also assume that the distribution of ^ is continuous, 

max{Ffc(Qfe-(-5)} < 1 and A„ < max j '^^^ f^J''' -(s + rh)}. Then: 
k k I rk[Qk) ) 

1. Any optimal allocation X* has a positive market order quantity M* > if ^ holds. 

2. Any optimal allocation X* has a positive limit order quantity L* >0 if ^ holds. 

3. // ([7])-([8]) hold for all exchanges j = 1, . . . , K , a necessary and sufficient condition for opti- 
mality of an order allocation X* is that it solves equations ([9]) -(10). 
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Proof: Proposition [2] implies the existence of an optimal order allocation X* £ C. First, we define 

Xm = (3,0, • • • ,0) and prove that X* ^ Xm by contradiction. If Xm were optimal in the problem 

Q it would also be optimal in the same problem with a constraint Lk = 0,k j , for any one j. In 

other words, the solution (S, 0) would be optimal for any one-exchange problem, defined by using 

2s + / + 

only exchange j. But by our assumption, there existsJ such that A„ < — (s + rj) and 

. . . ^j[Qj) 

Proposition ^ implies that (5,0) is not optimal for the J-th single-exchange subproblem, leading 

to a contradiction. 

The function v{X,^) is bounded for X £C and for all ^, differentiable with respect to M and 
Lk,k = 1, . . . , K for X € C\ {Xm} for almost all ^. Applying the same theorem as in the proof of 
Proposition [s] we conclude that V{X) is differentiable for X G C\ {Xm} and we can compute all 
of its partial derivatives by interchanging the order of differentiation and integration. The KKT 
conditions for problem Q and X G C\ {Xm} are 

s + f- \^¥{A{X*,0 <S) + XoHMX^O >S)-fi„ = (17) 
-(s + r,)P(e, > Qk + LI) - X^F{A{X\0 <S,^k>Qk + Ll) + 

\,nA{X\i)>S,ik>Qk + Ll)-iJik = ^, k = l,...,K (18) 

M>0, Lfc>0, ^o>0, Mfe>0, HoM = 0, /ifcLfe = 0, k = l,...,K (19) 



Since the objective function V{.) is convex, conditions (17)-(19) are both necessary and sufficient 



for optimality. The first result of this proposition follows from considering any X with M = 0: 



V{X) > X^SF [n{^k < Qk}j -{s + max{r,})5P (^Ql^ < Qk}j > (s + f)S = ^(Xm) and we 
already argued that 3X* with V{X*) < V{Xm), so X* / X and therefore M* > 0. Rearranging 



terms in a j'-th equality (18) we obtain 



> Q, + L*) [X, -{s + r,)- (A„ + X,)¥{A{X\^) < S\^, > Q, + L*)] - /i, = (20) 



The term in square brackets in (20) is negative for any X € C\ {Xm} with Lj = 0, because 



p(^(x,o < s\^, > Q, + L,) > p n {6 < Qk} 

by assumption and since fij > the condition (18) cannot be satisfied with L* = 0. We showed 



that M* > 0, L* > for all j = 1, . . . , K and therefore, fio = fix = ■ ■ ■ = fix = by complimentary 
slackness. Then the KKT conditions ([l7])-([l9]) reduce to ([9j)-([l0j). ■ 
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Corollary Consider the case of two exchanges with ^1,^2 that are independent and have continuous 
distributions. If 

1. max{Fk{Qk + S)}<l 

ro I # I ^ 2s + / + max{rfe} 

2. A„ < max <{ ^ , [s + rk) >, A„ > (s + max|rfc|) 



k=l,2 



FkiQk) 

3. Fi(Q0<l-^,i^2(Q2)<l '^'^ 



i^l(Ql)i^2(g2 



fe=l,2 



i/ien i/iere exists an optimal order allocation X* = (Af*,L*,L2) G mt{C} and it verifies (11a 



11c) 



Proof: Solutions on the boundary of C are sub-optimal: M* = and M* = S are ruled out by 
assumption 2, L\ = S — M and L2 = S — M are ruled out 



3y assumption 3 and (18). Solutions 



K 

with M* + Ll, = S are ruled out by directly checking 
ruled out by ( |18[ ). For example if L* = 0, then by Proposition [ 
complimentary slackness, F{A{X\0 < 5',C2 > Q2 + ^2) = > 5',6 > Q2 + ^2) = 0- But 

2>Q2+L*)>0- 



18). Finally, L* = and Ll = are also 
[i]M* + L* = 5' and in ([is]) ;U2 = by 



then (18) cannot hold because 



For any X G int{C}, A{X,S^) > 5" if and only if all the following three inequalities are satisfied: 



ei>Qi + S-M-L2 
6>Q2 + S-M-Li 
6+6>Qi + Q2 + 5-M 



(21a) 
(21b) 
(21c) 



These inequalities give a simple characterization of the event {A{X,^) > S} and their equivalence 
is directly verified by considering subsets of (^11^2) forming a complete partition of M^. 

Case 1: 6 > Qi + ^1,6 ><52 + ^2- Since Li + L2 + M > 5, we have A{X,^) = Li + L2 + M > S 



and at the same time all of the inequalities (21a 21c) are satisfied, so they are trivially equivalent 
in this case. 



Case 2: Qi + Li,Q2 < (,2 < Q2 + Because of the condition ^1 > Qi + Li, (21a) is satisfied. 



We have in this case that A{X,£,) = Li + ^2 - Q2 + M and thus A{X, ^) > 5 if and only if ( 21b ) is 



satisfied. Finally, ^1 > Qi + Li together with (21b) imply (21c), so A{X,£^) > S and (21a 21c) are 
equivalent in this case. 

Case 3: ^^2 > Q2 + -^2,Qi < < Qi + Li. Similarly to Case 2 we can show that inequalities 



(21a 21c) are satisfied if and only if A{X,^) > S. 

Case 4: Qi + S - M ~ L2 < S,i < Qi + Li,Q2 + S - M - < ^2 < Q2 + ^2- This set is non- 



empty because Q < S — M — Li < L2 and similarly for Li,L2 reversed. Inequalities (21a)-(21b) hold 
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trivially, only (21c) needs to be checked. We can write A{X,^) = — Qi + — Q2 + M > S ii and 



only if (21c) holds, so A{X,^) > 5 is equivalent to (21a 21c). 



Case 5: Outside of Cases 1-4, either (21a) or (21b) is not satisfied. If < Qi + 5* — M — L2,^,2 < 



Q2 + L2, then A{X,C) < S - M - L2 + L2 + M = S. The case 6 < Q2 + 5 - M - Li, ^1 < Qi + Li is 
completely symmetric, and it shows that neither A{X,^) > S nor ( 21a|21c ) hold in this case. 
Next, we use inequalities ( 21a|21c ) to characterize the set {A{X, ^) > S} in the first-order conditions 
([9])-(10). We observe that in the two-exchange case 



{AiX, > 5, 6 > Qi + Li} = {6 > Qi + Li, 6 > Q2 + 5 - M - Li} 

{A{X, >S,^2>Q2+L2}={C2>Q2 + L2,Cl>Ql+S-M- L2}, 

and then use the independence of ^1 and ^2 to compute 

F{A{X, > S|6 > Qi + L,) = F2{Q2 + S - M - L,) 
F{A{X, 0>S\i2>Q2 + L2) = h{Qi + S-M-L2) 

Together with ^ and (10), this leads to a pair of equations for limit orders sizes: 

A(Q2 + 5-M-L,) = V't"^' UQi + S-M-L2) = \^'X"' 

whose solution is given by L^^L^ from ( lla|llb" ). To obtain the equation (11c), we rewrite the first 
equation in ( 9|10 ) using the inequalities (21a 21c). Then P{A{X,^) > S) may be computed as the 
integral of the product measure Fi (8) F2 over the region defined by 



U{Q,S,M,Li,L2)={{Xi,X2)eR\ Xi> Q1+S-M-L2, X2> Q2+S-M-L1, Xi+X2> Q1+Q2+S-M}. 

This integral is given by 



P{A{X, > 5) = Fi ^ F2 {U{Q, S, M, Li, L2)) 

Ql+Li 

= F,{Q, + L,)F2iQ2 + S-M-L^) + 



MQi + Q2 + S-M- Xi)dFi(xi) 



K-{s + f) 



+S-A/-L2 
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